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The conversion of nuclear matter to quark matter in the 
core of a rotating neutron star alters its moment of inertia. 
Hence the epoch over which conversion takes place will be sig- 
naled in the spin-down characteristics of pulsars. We find that 
an observable called the braking index should be easily mea- 
surable during the transition epoch and can have a value far 
removed (by orders of magnitude) from the canonical value of 
three expected for magnetic dipole radiation, and may have 
either sign. The duration of the transition epoch is governed 
by the slow loss of angular momentum to radiation and is 
further prolonged by the reduction in the moment of inertia 
caused by the phase change which can even introduce an era 
of spin-up. We estimate that about one in a hundred pulsars 
may be passing through this phase. The phenomenon is anal- 
ogous to "bachbending" observed in the moment of inertia of 
rotating nuclei observed in the 1970's, which also signaled a 
change in internal structure with changing spin. 

97.60Lf, 97.60.Gb, 97.10.Cv 

The deconfined phase of hadronic matter called quark 
matter is believed to have pervaded the early universe 
and may reside as a permanent component of neutron 
stars in their dense high-pressure cores [0-|6) ■ However no 
means of detecting its presence has been found because 
the properties of neutron stars and those with a quark 
matter core are expected to be very similar. Alternately, 
instead of looking to the properties of the star itself, we 
study the spin-down behavior of a rotating star with the 
realization that changes in internal structure as the star 
spins down will be reflected in the moment of inertia and 
hence in the deceleration. 

Pulsars are born with an enormous store of angu- 
lar momentum and rotational energy which they radi- 
ate slowly over millions of years by the weak processes of 
electromagnetic radiation and a wind of electron-positron 
pairs [f7| — p~0|] . When rotating rapidly, a pulsar is centrifu- 
gally flattened. The interior density will rise with de- 
creasing angular velocity and may attain the critical den- 
sity for a phase transition. First at the center and then 
in an expanding region, matter will be converted from 



the relatively incompressible nuclear matter phase to the 
highly compressible quark matter phase. The weight of 
the overlaying layers of nuclear matter will compress the 
quark matter core and the entire star will shrink. The 
mass concentration will be further enhanced by the in- 
creasing gravitational attraction of the core on the over- 
laying nuclear matter. The moment of inertia thus de- 
creases anomalously with decreasing angular velocity as 
the new phase slowly engulfs a growing fraction of the 
star. 
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FIG. 1. Evolution of radial location between boundaries 
of various phases in a neutron star of given baryon number as 
a function of moment of inertia, which is a decreasing function 
of time. Stellar mass is 1.42M© for the nonrotating counter- 
part. Formation of a pure quark matter core occurs below 
/ ~ 107 km 3 (fl ~ 1400 rad/s). Note the decrease of radius 
with time. See Ref. |nj for a description of spatial geometry 
in mixed phase. 

The decrease of moment of inertia caused by the phase 
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transition is superposed on the natural response of the 
stellar shape to a decreasing centrifugal force occasioned 
by radiation loss. Therefore, to conserve angular mo- 
mentum not carried off by radiation, the deceleration Q 
of the angular velocity must respond by decreasing in 
absolute magnitude and may actually change sign. The 
pulsar may spin up for a time, just as an ice skater spins 
up upon contraction of the arms before air resistance and 
friction reestablish spin-down. 

Such an anomalous decrease of the moment of inertia 
as described is analogous to the phenomenon of "back- 
bending" in the rotational bands of nuclei predicted by 
Mottelson and Valatin and observed years ago [[l^-[l4||. 
The moment of inertia of a nucleus changes anomalously 
because of a change in phase from a nucleon spin-aligned 
phase at high angular momentum to a pair-correlated su- 
perfluid phase at low. The connection between the mo- 
ment of inertia and the internal structure of a neutron 
star is shown in Fig. [L and between moment of inertia 
and frequency in Fig. 
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FIG. 2. Moment of inertia as a function of rotational an- 
gular velocity over the epoch during which an increasing cen- 
tral volume of the star is passing into the deconfined quark 
matter phase. The temporal development is from large to 
small /. Spin-up of the pulsar is clearly evident. An analo- 
gous phenomenon known as "backbending" is observed in the 
rotational properties of nuclei. 



The property of asymptotic freedom of quarks makes 
a phase change from confined to deconfined matter in- 
evitable for sufficiently high energy density of hadronic 
matter. We have therefore referred to the transition as 
being from the confined to deconfined phase. However, 
it is apparent in our discussion and results, that the phe- 
nomenon we discuss is not unique to this particular tran- 
sition. It is simply the most plausible in our view, but 
any phase transition in dense hadronic matter is inter- 
esting. 



Now we develop the appropriate measure for detecting 
the occurrence of a phase transition and its envelopment 
of an increasing proportion of the mass and volume of 
the star with time. The production of electromagnetic 
radiation and a wind of relativistic pairs by the rotating 
magnetic field of the star exert a torque on the star which 
is characteristic of magnetic dipole radiation. The energy 
loss equation representing processes of multipolarity n is 
of the form 
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where n = 3 for magnetic dipole radiation. The rate of 
change of frequency is governed by 
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For low frequency or if changes in / are ignored this re- 
duces to the usual form quoted in the literature. O = 
-K O n where K = C/I, (cf. Ref. HQ). 

The dimensionless measurable quantity OO/O 2 is re- 
ferred to as the "braking index". It would be equal to 
the intrinsic index n of the energy-loss mechanism (|l|) 
if the frequency were small or if the moment of inertia 
were a constant. However these conditions are not usu- 
ally fulfilled and the measurable quantity is not constant. 
Rather it has the value 
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where I' = dl /dO and I" = dl 2 /dfl 2 . The progression of 
the new phase through the central region of the star will 
be signaled by an anomalous value of the braking index, 
far removed from the canonical value of n. 

Because the pulsar rotational energy is coupled to 
weak processes, O is small and none of the quantities in 
([|) will change appreciably over any observational time 
span. However the signal is carried in non-zero deriva- 
tives dl/dVt and dl 2 /d£l 2 ; these are large during the 
phase transition epoch because of the progressive con- 
version of nuclear matter into compressible quark matter. 
As can be seen from (^) , large derivatives of the moment 
of inertia will produce enormous deviations of the brak- 
ing index from its canonical value while the region of 
the star occupied by quark matter is growing. Since the 
growth is paced by the slow spin-down of the pulsar, we 
will find that the signal is "on" over an extended epoch. 

The behavior of the moment of inertia in the critical 
frequency interval for our model star (which is described 
later) is shown in Fig. ^. As the pulsar evolves in time 
(decreasing I) the derivative dl/dfl passes through two 
singularities, switching between +oo and -co at each 
turning point of O. From (|^) it is clear that the deceler- 
ation O will pass through zero and change sign at both 
turning points becoming an acceleration in the central 
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part of the 'S'; the pulsar spins up for a time. Moreover, 
—I" has similar singularities as can be found from 

\dnj di 2 [ ' 

Consequently n(O) goes to ±00 at the two turning points 
respectively as shown in Fig. ||. We have plotted the 
braking index as a function of / because / decreases 
monotonically with the time, unlike f2. 

(The spin-up referred to above has nothing to do with 
the minuscule spin-up known as a pulsar glitch. The rel- 
ative change in moment of inertia in a glitch episode is 
very small (A/// ~ — Af2/f2 ~ 10~ 6 or smaller) and ap- 
proximates closely a continuous response of the star to 
changing frequency on any time scale that is large com- 
pared to the glitch and recovery interval. Excursions of 
such a magnitude as quoted would fall within the thick- 
ness of the line in Fig. ||.) 

A change in moment of inertia owing to a change in 
phase such as we have described is evidently a robust 
phenomenon. However, we cannot be sure that nature 
will respond as strongly as our model does. Backbending, 
as in Fig. ^, is an extreme response of the moment of 
inertia to the progression of a phase transition through 
the central region of the star. Instead, the transition of 
the moment of inertia from that of a neutron star (at 
high f2) to a hybrid star may be a single- valued function 
of Q. However the transition between these two moments 
can still be marked by a large first derivative /' and a 
second derivative I" that is large and changes sign in the 
transition interval. In this case, the braking index will 
not swing between ±00 but nonetheless can attain large 
positive and negative values. 

Typically it is difficult to measure Q (and hence the 
braking index) because of timing noise. Only four brak- 
ing indices are presently known. However, for a star that 
is passing through the phase transition epoch, the decel- 
eration (||) is reduced markedly (even changing to accel- 
eration). The second derivative must therefore be large 
in absolute magnitude through all but the central por- 
tion of the epoch where spin-up occurs. Consequently, 
the second derivative of frequency should be easier to 
measure for pulsars passing through the epoch than for 
typical pulsars. Hence, difficulty in measuring the second 
derivative may be used as a de-selection criterion. 

We estimate the plausibility of observing in the pul- 
sar population a signal of the kind that we find in our 
calculation. Assume that neutron stars are created in a 
narrow mass window (as present evidence suggests). The 
duration over which the observable index is anomalous is 
AT w — A57/57 where Af2 is the frequency interval of the 
anomaly. The range in which n is smaller than zero or 
larger than six (Fig. ||) is Ail w 100 rad/s. Take a typi- 
cal period derivative of P ~ 10~ 16 s/s to find AT ~ 10 5 
years. The signal would endure for 1/100 of a typical 



active pulsar lifetime. Similarly we estimate that spin- 
up would last for about 1/6'th of that time. Given that 
more than 700 pulsars are presently known (200 of which 
have been discovered in the last several years) about 7 of 
these may be signaling the growth of a central region of 
new phase. 
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FIG. 3. The braking index as a function of moment of 
inertia signals by its large departures from n = 3 the pro- 
gression of the deconfined phase in the stellar interior. The 
braking index would be identically 3 if it were not for the 
response of / to rotation, and in any case tends to 3 for low 
frequency (for magnetic dipole radiation). 

We describe briefly the calculation. The equations de- 
scribing the configurations of rotating stellar structures 
are solved for a sequence having the same baryon number 
but different rotational frequencies. The usual expres- 
sion for the moment of inertia in General Relativity due 
to Hartle is not adequate for our purpose (see Ref. JIt]] 
equations (49) and (62)). Not only are effects of internal 
constitutional changes in the star absent but also absent 
are the effects of the alteration of the metric of space- 
time by rotation, the dragging of local inertial frames, 
and even centrifugal flattening. Instead, to calculate the 
moment of inertia, we must use an expression that incor- 
porates the above effects as derived by Glendenning and 
Weber @,||. 

The stellar model is described as follows. Neutron star 
matter has a charge neutral composition of hadrons con- 
sisting of members of the baryon octet together with lep- 
tons when in the purely confined phase. The proper- 
ties of such matter are calculated in a covariant mean 
field theory as described in Refs. pp|-|2^]. The values 
of the parameters that define the coupling constants of 
the theory are certain fairly well constrained properties 
of nuclear matter and hypernuclei as described in the 
references; (binding energy of symmetric nuclear matter 
B/A = -16.3 MeV, saturation density p = 0.153 far 3 , 
compression modulus K — 300 MeV, symmetry energy 



3 



32.5 MeV, nucleon effective mass at 
0.7m and ratio of hyperon to nucleon 



that yield, together 



coefficient a sym 
saturation m* at 

couplings x a = 0.6, — 0.653 
with the foregoing parameters, the correct A binding in 
nuclear matter pH). Quark matter is treated in a ver- 
sion of the MIT bag model with the three light flavor 
quarks (m u = = 0, m s = 150 MeV) as described in 
Rcf. |||. A value of the bag constant B 1 ' 4 = 180 MeV 
is employed. (See Refs. [^2 26 27]] for a correct treatment 
of first order phase transitions in multi-component sub- 
stances such as neutron star matter for which baryon and 
electric charge are the independent conserved charges.) 
Very little is known about the high-density properties of 
matter and our calculation does not imply a prediction of 
the rotational frequency or stellar mass at which a phase 
transition will occur. Rather it shows what the signal 
could be if conditions are attained for the phase change. 

So far as we know, there is little difference in the prop- 
erties of a neutron star that has no quark core and one 
that has already fully developed one. But a strong signal 
may be associated with the gradual conversion of matter 
from one phase into the other as the conversion is paced 
by the slow loss of angular momentum in the processes 
of electromagnetic radiation and electron-positron wind. 
The important observational features by which such an 
epoch could be identified are: 

(1) The braking index has a value far from the canonical 
value, possibly by orders of magnitude and can be of ei- 
ther sign. 

(2) The epoch over which the braking index is anomalous 
is long because pulsars spin down slowly under ordinary 
circumstances but even more slowly when |/'| is large (see 
(§)■ 

(3) The pulsar may be observed to be spinning up. (To 
avoid confusion with spin-up due to accretion, only iso- 
lated pulsars are relevant.) 

(4) Except for the central part of the spin-up era, the 
derivative fl is large and therefore easy to measure and 
so also is the braking index. 

(5) Difficulty in measuring ft can be used to deselect 
phase transition candidates. 

(6) An estimated 1/100 pulsars may be passing through 
the transition epoch. 

Pulsar observations are still in their infancy. It takes a 
considerable time-span of data to measure the braking in- 
dex. And many of the presently known pulsars have only 
recently been discovered. We conclude from our work 
that it is plausible that the phase transition signal can 
be observed. It would be a momentous discovery to find 
that a phase of matter that existed in the early universe 
inhabits the cores of some neutron stars. 
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